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NOTATION 


2 

Aj^  “ Recoil  rod  area,  in 

2 

Aq  “ Control  rod  area,  in 

2 

Ap  - Floating  piston  area.  in 

2 

Ajg  - Total  area  of  floating  piston  = A^  + A^.  in 

Aj^  - Area  of  orifice  between  recoil  and  recuperating 

2 

chamber,  in 

A2  - Area  of  orifice  between  recuperating  and  control 
chamber,  in^ 

A^(x)  - Variable  area  of  the  groove  in  the  floating 

2 

piston  at  position  x.  in 
B(t)  - Breech  force  at  time  t.  Ibf 
C-^  - Discharge  coefficient  for  orifice  A^^. 

C2  “ Discharge  coefficient  for  orifice  A2. 

~ Discharge  coefficient  for  orifice  A^. 

Cg  - Equivalent  friction  coefficient  for  frictional 
loss  at  orifices 

Fj^  - Dry  friction  at  recoil  piston. 

Fp  - Dry  friction  at  floating  piston. 

F - Equivalent  dry  friction, 
g^^  - Position  feedback  gain. 
g2  - Velocity  feedback  gain. 

g^^g^  ” other  feedback  gains  for  nonlinear  feedback. 

J - objective  function 


- Components  of  objective  function 


2 

mp  - Mass  of  the  floating  piston  Ibf  sec  /in 

2 

mj^  - Mass  of  the  recoiling  parts  Ibf  sec  /in 

2 

m - Equivalent  mass  Ibf  sec  /in 

q 

2 

Pq  - Initial  gas  pressure  Ibf/in 

2 

Pq  - Gas  pressure  at  time  t Ibf/in 
R - Gas  constant 

RDPL(t)  - Rod  pull  at  time  t Ibf 

RDPLD  (t)  - Desired  Rod  pull  at  time  t Ibf 

T - Recoil  time  sec 

2 

U(t)  - Open  area  of  servo  valve  in 
X - Position  of  recoiling  parts  with  respect  to 
recoil  rod  in 


X - Velocity  of  recoiling  parts  with  respect  to 
recoil  rod  in/sec 
X - Acceleration  of  recoiling  parts 
y - Position  of  floating  piston  with  respect  to 
recoiling  parts  ^ 

w - Total  weight  of  the  recoiling  parts. 

Wg  - Equivalent  weight  of  recoiling  parts. 


X 


w. 


w. 


w. 


w. 


- Weighting  factors  for  components  in  the 
objective  function. 


a - The  angle  of  elevation  of  the  gun  barrel 
p - Mass  density  of  the  hydraulic  fluid. 


radians 
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SUMMARY 


Mathematical  model  for  a conventional  hydro- 
pneumatic recoil  mechanism  is  developed  from  physical 
laws.  This  mathematical  rtiodel  is  simulated  on  a digital 
computer  by  Phase-Plane-Delta  method.  The  time  histories 
of  all  pertinent  parameters  such  as  position  of  recoiling 
parts,  hydraulic  pressures  in  different  chambers,  rod 
pull  are  available  and  can  be  plotted. 

A linear  state  feedback  control  system  is  proposed 
to  adapt  this  conventional  recoil  mechanism  to  perform 
satisfactorily  for  all  firing  zones.  The  control  and 
optimi-zation  problem  is  formulated.  An  objective  function 
with  direct  relations  to  the  performance  and  constraints 
of  the  problem  is  formulated.  Davidon-Fletcher-Powell 
variable  metric  method  of  nonlinear  optimization  is  used 
to  design  the  feedback  gains  of  the  state  feedback  con- 
trol law. 

The  above  optimisation  procedure  is  applied  to  M-  37 
recoil  mechanism  for  zone  5 through  8.  This  recoil 
mechanism  originally  was  designed  for  zone  7.  One  set 
of  feedback  gains  is  designed  for  each  zone  and  switching 
of  control  laws  is  suggested.  A significant  improvement 
in  the  recoil  force  time  trajectory  shape  and  reduction 
in  peak  force  of  2.5  to  25  percent  is  obtained. 
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INTRODUCTION 

Recoil  mechanisms  dissipate  energy  of  the  reaction 
of  gunfire  at  a controlled  rate  so  as  to  minimize  the 
recoil  force  transferred  to  the  carriage  of  a weapon 
system  without  exceeding  available  recoil  length.  Nerdahl 
and  Frantz  t 9 ] have  developed  three  degrees  of  freedom 
nonlinear  models  of  hydropneumatic  recoil  mechanisms  and 
a procedure  to  design  a variable  area  orifice  to  control 
the  energy  dissipation.  The  procedure  defines  a control 
function  or  desired  control  recoil  force-time  trajectory 
for  a design  firing  charge  and  with  the  help  of  digital 
simulation  of  the  model, computes  the  orifice  area  at 
different  positions  of  the  recoil  mechanism.  This 
design  performs  satisfactorily  for  the  designated  firing 
charge,  but  far  fi^om  optimum  for  other  firing  charges. 
Thus,  a control  system  which. can  adapt  to  different  firi^ig 
aones  is  desirable. 

A linear  state  feedback  control  system  with 
variable  gains  is  proposed  in  the  report.  A separate 
control  law  is  designed  for  each  firing  charge  and  the 
control  law  corresponding  to  the  charge  being  fired  is 
selected  from  this  predesigned  set.  This  control  scheme 
can  be  implemented  by  adding  a servo  valve  operating  in 
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tandem  with  the  variable  area  orifice.  The  feedback 
gains  for  the  servo  valve  can  be  selected  from  a 
predesigned  set  by  identifying  the  charge  being  fired 
by  sensing  signals  such  as  acceleration  with  the  help 
of  a microprocessor  or  special  purpose  digital  electronics. 

The  design  of  the  optimal  control  law  is  complicated 
by  the  highly  nonlinear  nature  of  the  model  induced  by 
turbulent  flow  through  orifices,  dry  friction  at  piston- 
cylinder  surfaces  and  adiabatic  gas  compression.  The 
nonlinear  second  order  model  is  simulated  on  a digital 
computer  using  Phase-plane-delta  method.  A multifactor 
objective  function  with  direct  physical  interpretation 
is  developed  as  a function  of  the  variables  computed  in 
the  simulation.  Davidon-Fletcher-Powell  variable  metric 
quasi-Newton  nonlinear  optimisation  algorithm  is  used 
to  minimize  the  objective  function.  The  procedure  is 
applied  to  a M-37.1Q5  mm  recoil  mechanism.  Also,  velocity 
feedback  for  lower  zones  1,  5,  and  6 is  evaluated. 
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I 

MATHEMATICAL  MODEL  FOR  HYDROPNEUMATIC  RECOIL  MECHANISM 

A schematic  diagram  for  a conventional  hydropneu- 
matic recoil  mechanism  is  shown  in  Fig.  1.  A mathematical 
model  based  on  the  physical  laws  is  developed  with  the 
following  assumptions. 

1)  Flow  through  orifices  is  a potential  flow. 

2)  Temperature  variations  do  not  affect  the  discharge 
coefficients  for  the  orifice. 

3)  The  recoil  mechanism  is  secured  to  the  carriage 
through  a rigid  link  called  the  rod  and  the 
carriage  rests  on  a rigid  support. 

4)  Only  translation  of  the  recoil  parts  in  direction  of 
firing  is  considered. 

5)  There  is  no  cavitation  in  any  cheimber. 


PRESSURE  RELATIONS 


« T,  /P\  / R \2  .2 

Pr  - ^ 


2 

= i sgn  X 

D ^R  2 

where  T^  = (f)  (^) 


Pi  - P2  = (7-S-) 


C .2  .2 


2'  'A^C-  ^ 


,p,  , •^C^R  ,2  .2 

= (7)  (a  r A ^ ^ ^ 

^2''2% 


where  T, 


_ • 2 • 

= T2  X sgn  X 

<2>  'AjC^A^’ 


'1  • '’3  = 'f> 

' <f> 


= T_  X sgn  X 


where  T3  = (f)  ^ 


(1) 


(2) 


(3) 
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The  gas  pressure  relation 


V, 


p = p {— 

G 0 ‘V„  - Aj^x 


0 ,R 


(4) 


Thus  the  pressure  equations  are 


Pj^  = + T^  X sgn  X 


1 " "1 


.2 


P2  = Pj^  - T2  X sgn  X 


P^  = Pj^  - Tj  X sgn  X 

P = p ( 5 

G 0 Wq  - Aj,x^ 


(1) 

(2) 

(3) 

(4) 


Force  balance  for  floating  piston 


sgn  X 


(5) 


Substituting  for  Pj^,  P2»  P3  and  P^  from  equations  1,  2 
3,  & 4 in  5 we  have 


* 2 * 

mp(l+^)  X = mpg  sin  a + (P^^-T^x'^sgnx) 


2 ^0 

+ Aj,{P^-T2i  sqni)  - A^Pp  sgni 


= mg  sin  a - (ApT2+A^T2)  x sgnx 


V 


0 vR 


+ (Aj^+A^)Pi  + m g sin  a - A^Pq  -)  - F sgnx 


0 "R" 
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R.  .. 


Aj^Pl  = mp(l+^)  X - iHpg  sin  a - (ApT3+A^T2) sgnx 

(6) 


* ^i'’o  ‘ V^’  ^ 


Force  balance  for  recoil  mass 


nij^x  = B(t)  + mj^g  sin  a - A^P^  - 


+ FpSgn^  - Aj^Pj^  - Fj^sgni 


Adding  Equations  (5)  and  (7)  we  have 


(7) 


[m_  + m (1+^)  ] X = (m  +mp)g  sin  a + B(t)  - 

K p Aj^  p K 

Vr  ■ ''r=9">' 


(8) 


Substituting  for  P^^  in  Eq.  (8)  from  Eq.  (6)  we  have 


R 


[n'R+ro  (1+— ) ] X = (m  + mj^)g  sin  a + B(t) 


'N 


“Fj^sgnx  - Aj^Tj^x  sgnx  - — [mp(l+^)x  - rrVpg  sina 


V 


+ (A^T2  + A15T3)  x^sgnx  + A^^Pq  ^ + FpSgnx] 


[m  +m  (1+— ) ] g sin  a + B(t) 
« P An 


«p  o 

[ArTi+^(AcT2+ApT3)  ] x"=sgnx 


R. 


- A P {— ^ 


Ap  *• 
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Finally 

[mg+mp(l+j|)^lX  + =9"i 


V 


0 "R 


N 


= [m„+m  (H~)  ] g sin  a + B(t) 


R p A, 


N 


(9) 


Let 


niq  - equivalent  mass 
= "r  + % (1  + 


C - equivalent  damping  coeffienct 

= ' Vl  ^ ^ <'^c'f2  * Vs'  > 


Fg  - equivalent  dry  friction 
^^R  ^p 

W - equivalent  weight 

Hi 

Ar 

= Imj^  + mp  (1  + ^)  ] g sin  a 
Thus  equation  (9)  can  be  rewritten  as 


mqX  + CqX  sgnx  + FqSgnx  + Aj^Pq(^  - JlR 

V„ 


= Wq  + B(t) 


(10) 


Note  - C is  not  a constant  and  is  a function  of  the 

q 

position  of  the  floating  piston  due  to  the  variable 

3* 


orifice  area  A 
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The  rod  pull  is 

RDPL  = Pj^Aj^  + Fj^sgnx  (11) 

From  equations  (8)  and  (11)  we  have 

RDPL  = ^ a + B(t) 

Thus  equation  (10)  is  the  final  model  of  the  form 

mx  + (a+f  (x)  )x^sgnx  + k(l-aQX)”^=  W + B(t) 
and  equations  1,  2,  3,  4,  and  12  are  relations  for  all 
pertinent  variables. 


(12) 


10 


II 

DIGITAL  SIMULATION  OF  RECOIL  MECHANISM  MODEL 

The  nonlinear  second  order  model  developed  in 
Section  I Equation  (10)  can  be  simulated  on  a high  speed 
digital  computer  by  phase-plane-delta  method.  This 
method  transforms  a forced  nonlinear  model  into  a linear 
oscillator  for  a short  interval  of  time  and  is  explained 
in  detail  in  Appendix  I.  The  simulation  results  in 
time  histories  of  all  the  pertinent  variables  such  as 
position,  velocity,  rod  pull,  and  all  the  pressures. 

This  method  is  computationally  very  efficient. 

The  model  in  equation  (10)  can  be  reformulated  in 
phase-plane-delta  format  as  follows 

m X + c^x^sgni  + A^P^  (1-^x)  ^ + F^sgnS:  = W +B(t) 

4 4 ^ ^ ''o  ^ ” 

or 

X + p^  (x  + 6x)  = 0 (13) 

where 

6x  = -X  + -~lC  x^sgn*  + A P-(l-^x)  ^ 

P m ^ 0 

^ (14) 

+ F sgnx  - - B(t)} 

4 q 

From  initial  conditions  Xq  = 0 , x^  = 0 , the  simulation 
is  started  and  run  until  the  recoil  mechanism  starts 
counter-recoil. 

A fortran  program  is  written  to  simulate  the  recoil 
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mechanism  based  on  equations  (13)  and  (14) . 

An  example  simulation  was  run  for  M-37  recoil 
mechanism  zone  7 firing  charge.  The  phase  plane, 

X,  X,  Pq,  and  rod  pull  are  plotted  in  Figures  2,  3, 

4,  5,  6,  and  7. 

The  simulation  results  are  in  agreement  with  a 
digital  simulation  run  by  the  Ware  Simulation  Division 
at  Rock  Island  Arsenal  using  Continuous  System  Modeling 
Program  (CSMP) . 


1 1 1 1 1 1 1 1 mu  H i i l ij  I i 1 1 1 1 1 1 1 1 1 1 1 1 1 I I 1 1 1 1 1 1 I t 1 1 1 1 i 1 1 1 1 1 i H 
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Figure  2 Phase  plane  for  Zone 
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Figure  5 Cavitation  Pressure  for  Zone 
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Figure  6 Gas  Pressure  for  Zone 
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Figure  7 Rod  Pull  for 
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III 

FORMULATION  OF  ADAPTIVE  CONTROL  AND 
OPTIMIZATION  PROBLEM 

The  control  groove  machined  in  the  floating  piston 
is  designed  for  a designated  zone  and  performs  satis- 
factorily when  used  for  the  design  zone.  The  purpose 
of  this  research  is  to  adapt  the  recoil  mechanism  such 
that  it  adapts  to  the  charges  being  fired  so  as  to 
perform  satisfactorily  without  violating  the  system 
constraints.  A modified  design  of  the  recoil  mechanism 
with  a servo  valve  operating  in  tandem  with  the  variable 
area  groove  is  shown  in  Figure  8.  The  area  of  the 
servo  valve  that  is  open  for  the  flow  of  hydraulic  fluid 
is  controlled  by  a feedback  law.  The  feedback  law  can 
be  changed  by  sensing  the  zone  being  fired.  This  adds 
the  flexibility  to  change  and  adapt  to  the  firing  zone. 

3.1  Linear  State  Feedback 

A linear  state  feedback  control  is  proposed  to 
control  the  area  of  the  servo  valve  and  is  of  the  form 
u = gj^x  + (15) 

where  g^^  and  g2  are  feedback  gains  and  u is  the  area 
of  the  servo  valve.  The  corresponding  modification  in 
the  mathematical  model  is  in  the  equivalent  damping 
coefficient  as  follows 


C 


and 


(16) 


The  modified  form  of  the  model  is 


• • 2 • *“R 

mx  + [a+f (x) +u (x,x) ] x sgnx  + k(l-aQx)  + Fsgni 


= W + B(t) 

The  task  of  finding  optimal  values  and  g2  of  the  feed- 
back gains  is  formulated  as  follows. 

Find  optimal  values  of  g^^and  g2  such  that  an 
objective  function  J(gj^,g2)  is  minimized  subject  to 
the  following  constraints.  The  system  follows  the  model 

1)  mgX+CgX^sgnX+FgSgnx+Aj^PQ(l-~  x)"^ 


where 


A„T,  + 


2 


2)  X 


max 


the  maximum  recoil  length 


£ Xj^,  the  available  recoil  length 


3)  There  is  no  cavitation  in  chamber  3,  i.e., 


P3(t)  > P 


3 min 


4)  The  maximum  servovalve  area  u = g^^x  + g2X 

U < U 

max  — m 
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3.2  Development  of  an  Objective  Function 

The  basic  performance  criterion  is  that  the  actual 
rod  pull  trajectory  should  follow  as  closely  as  possible 
the  desired  control  trajectory  for  every  zone.  So  let 
the  control  trajectory  be  RDPLD(t) 

Then , 


will  represent  the  integrated  error  or  least  square  error 
criterion.  To  be  able  to  use  unconstrained  optimization 
techniques,  penalty  functions  can  be  added  to  take  care 
of  constraints  and  hence  the  penalty  functions  are 


J(gi,g2)  = [RDPLD(t)  - RDPL(t)]^  dt 


T 


T 


T 


where 


W2  = 0 P3(t)  > P 


3min 


«3  = 0 U(t)  < 


X < X, 


m 


0 
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So  the  composite  objective  function  is 
T 

(RDPLD(t)  - RDPL(t))^  + W2(P3(t)  - P3  ^ 

+ w^(U(t)  - + w.(x-x„)^  } dt 

■i  m 4 m 

By  choosing  W2»  w^,  large  postive  numbers,  the  constraint 
violations  can  be  avoided. 

The  VAIOA  optimization  routine  availible  on  UNIVAC  1110 
at  UWMACC  which  uses  Davidon  - Fletcher  - Powell  (Appendix  III) 
method  was  used  to  find  optimal  values  of  g^^,  g2  for  zone  7 and 
results  are  shown  in  figure  9,  10.  The  control  trajectory 
was  chosen  arbitrarily.  Figure  9 is  with  no  feedback 
and  10  is  with  optimal  feedback  gains  of  g^^  = .2742E-3 
and  g2  = .354E-7.  The  actual  trajectory  follows  the 
desired  trajectory  in  Fig.  10  more  closely  than  in  Fig.  9. 

The  procedure  was  repeated  for  other  zones  with 
arbitrary  control  trajectory.  The  optimization  procedure 
is  very  sensitive  to  the  definition  of  the  control  tra- 
jectory and  thus  redefinition  of  the  objective  function 
without  the  control  trajectory  is  necessary.  The  main 
objective  can  be  reiterated  as  the  minimization  of  rod 
pull  and  a flat  trajectory  to  avoid  sudden  or  sharp 
changes  which  induce.-,  fatigue  failure.  The  first  deri- 
vative represents  the  smoothness  or  flatness  of  the 
trajectory  to  some  extent  and  hence,  a penalty  for  large 
first  derivatives  can  be  added  to  the  objective  fvinction. 


26000 
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Figure  9 Rod  Pull  for  Zone  7 with  IIo  Feedback 


Figure  10  Rod  Pull  for  Zone  7 With  Linear  Feedback 
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The  new  composite  objective  function  can  be  written  in 
five  parts  as 

T 

Jf  = (RDPL(t)-  C)^  dt 

where  C is  a constant  = 4200. 

Penalty  for  non-smoothness 
T 

^2  = K 'al 

Penalty  for  constraint  violation  for  cavitation 
J3  = (PjCt)  - P3  at 

Penalty  for  maximum  servo-valve  area  constraint 

J,  = l1  (U(t)  - at 

Penalty  for  maximum  recoil  length  constraint 

J5  ‘ 1]  « - 

The  composite  objective  function  being 


with  w^  = 

0 

if 

P^(t)  > P^ 

j J min 

II 

0 

if 

U(t)  < U„ 
— m 

II 

in 

0 

if 

X < 

m 

The  weighting  factors  w^^  and  W2  provide  a trade  off 
between  the  minimum  error  and  flatness  of  the  rod  pull 
trajectory.  Large  w^^  with  respect  to  W2  will  result  in 
sharp  trajectory  with  a peak  in  initial  stage  of  recoil. 
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Large  W2  with  respect  to  will  result  in  a sharp  peak  at 
a later  stage  in  recoil  and  also  possibly  will  result  in 
cavitation  at  variable  area  orifice  3.  By  choosing  w^^ 
and  W2  in  between  these  extremes/  a satisfactory  shape  of 
the  trajectory  can  be  achieved. 

The  weighting  factors  w^/  w^,  and  w^  emphasize  or 
de-emphasize  the  penalty  functions  for  constraint  viola- 
tions. If  any  constraint  is  violated,  the  corresponding 
weighting  factor  can  be  increased  so  as  to  force  the 
optimization  algorithm  to  choose  a feasible  solution. 
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3.3  Results  for  M-37  Recoil  Mechanism 

The  linear  state  feedback  control  system  and 
optimization  procedure  discussed  in  Sections  3.1  and 
3.2  was  applied  to  M-37  recoil  mechanism.  M-37  was 
designed  for  a firing  charge  zone  7.  The  design  data 
and  breech  forces  for  zones  1,  5,  6,  1,  and  8 are 
tabulated  in  Appendix  III  and  was  provided  by  the  Ware 
Simulation  Division  at  Rock  Island  Arsenal.  The  breech 
forces  used  are  from  simulated  breech  forces.  The  breech 
force  for  zone  8 was  arrived  at  by  multiplying  breech  force 
for  zone  7 by  1.2  due  to  unavailability  of  data. 

The  feedback  gains  for  zones  5,  6,  7,  and  8 and  other 
parameters  are  presented  in  Table  1.  The  rod  pull 
characteristics  are  plotted  in  Figures  11  through  20  with 
no  feedback  and  with  optimal  feedback. 

Figures  11  and  13  portray  the  time  history  of  the 
rod  pull  for  no  feedback  and  optimum  linear  feedback 
control  for  zone  8.  Zone  8 was  arrived  at  by  multiplying 
zone  7 breech  force  by  1.2.  Though  the  trajectory  for 
no  feedback  looks  very  satisfactory,  the  Figure  12 
which  plots  the  pressure  P3(t)  reveals  that  there  is 
cavitation  and  hence  the  model  is  no  longer  valid ► 

Figure  13  shows  trajectory  not  very  flat  but  there  is  no 
cavitation  as  shown  in  Fig.  14.  This  is  achieved  through 
optimization  with  penalty  for  cavitation.  With  Random 


Table  1 Results  of  Optimum  Feedback  for  Zones  5,  6/  7 and 
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Search  technique  was  increased  up  to  20  psi. 

With  reasonable  cavitation  can  completely  be 

eliminated. 

Figures  15  and  16  are  rod  pull  trajectories  for 
zone  7 with  no  feedback  and  with  optimal  feedback.  The 
improvement  is  not  significant.  The  improvement  was 
significant  as  shown  in  Figures  9 and  10  even  for  zone  7 
with  the  first  objective  function.  The  relative  iit^rove- 
ment  in  peak  force  is  about  2.5%  and  recoil  length  is 
longer  by  0.6  inches.  There  is  no  change  in  the  total 
recoil  time  of  .13  seconds. 

The  rod  pull  trajectory  for  zone  6 with  no  feedback 
is  a more  triangular  with  maximum  rod  pull  of  14,700 
lbs.  (Figure  17).  The  optimal  feedback  control  law  results 
in  a very  flat  trajectory  with  maximum  force  of  12,900 
which  is  a significant  reduction  of  about  12%  (Figure  18) . 
The  recoil  length  is  increased  from  25.8  inches  to 
27.27  inches  but  recoil  time  is  reduced  from  .17  seconds 
to  .158  seconds. 

The  rod  pull  trajectory  for  zone  5 with  no  feedback 
is  a sharp  triangular  one  with  peak  of  11,360  lbs.  (Figure 
19).  The  optimal  feedback  control  law  reduces  this  force 
to  8500  lbs  with  a flat  and  trapezoidal  trajectory.  The 
recoil  length  and  time  are  both  increased.  The  percentage 
reduction  in  recoil  force  is  25%. 


niiiiiliiiiiiiiiliiiijiijiljujiijiijiiniimiiimiiiri 
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Figure  15  Rod  Pull  for  Zone  7 With  No  Feedback 
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Figure  Ig  Rod  Pull  for  Zone  7 with  Feedback 
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Figure  17  Rod  Pull  for  Zone  6 with  No  Feedback 
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Figure  18  Rod  Pull  for  Zone  6 with  Feedback 
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3.4  Tachometer  Feedback 

To  investigate  different  control  strategies,  only 
velocity  feedback  of  the  form  u = g^^x  was  studied.  In 
general,  it  worked  very  well  for  lower  zones  6,  5 and  1 
but  failed  to  improve  the  trajectories  for  zone  7 and 
made  it  worse  for  zone  8.  The  optimum  feedback  gain 
was  arrived  at  by  trial  and  error. 

Trajectory  for  zone  7 with  tachometer  feedback 
gain  g^^  = .006E-3  (Fig.  21)  is  not  much  of  an  improve- 
ment over  one  without  feedback  (Fig. 15).  Fig.  22  for 
zone  6 with  g^^  = .025E-3  shows  a substantial  improvement 
over  Fig.  17,  the  maximum  force  being  13,200  lbs.; 
though  it  is  not. better  than  linear  state  feedback  (Fig. 18). 
Fig.  23  for  zone  5 with  g^^  = .07E-3  is  trapezoidal  and 
significant  reduction  in  maximum  rod  pull  to  9,400  lbs., 
but  is  not  better  than  linear  state  feedback  of  Fig.  20. 
Trajectory  for  zone  1 with  g^^  = .35E-3  (Fig.  24)  is  also 
very  good  with  reduction  to  5,350  lbs.  over  Fig.  25  (24%). 

• Non-linear  feedback  control  laws  of  the  form 

U = gj^x  + g^x  + g3X^  + g^k^ 

were  investigated  for  zone  7.  The  trajectory  obtained 
were  with  no  improvement  and  sometimes  were  worse  than 
those  with  no  feedback. 
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Figure  21  Rod  Pull  for  Zone  7 with  Velocity  Feedback 
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Figure  22  Rod  Pull  for  Zone  6 with  Velocity  Feedback 
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Figure  23  Rod  Pull  for  Zone  5 with  Velocity 
Feedback 
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Figure  24  Rod  Pull  for  Zone  1 with  Velocity 
Feedback 
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Figure  25  Rod  Pull  for  Zone  1 with  No  Feedback 
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CONCLUSIONS 

A mathematical  model  for  a conventional  hydropneumatic 
recoil  mechanism  was  developed  and  simulated  on  a digital 
computer  by  method  of  Phase-Plane-Delta.  A linear  state 
feedback  control  system  was  proposed  which  can  be  imple- 
mented by  retrofitting  the  present  designs  with  a servo- 
valve to  operate  in  tandem  with  the  variable  area  groove 
in  the  floating  piston. 

An  objective  function  with  direct  relation  to 
performance  and  physical  constraints  of  the  system  was 
developed.  Davidon-Fletcher-Powell  nonlinear  optimization 
algorithm,  was  chosen  to  optimize  this  objective  function. 
The  procedure  was  applied  to  M-37.  recoil  mechanism  with 
reduction  in  peak  recoil  forces  from  25  to  2.5%  for  lower 
zones  and  cavitation  was  avoided  for  zone  8.  Tachometer 
feedback  was  shown  to  be  effective  for  low  zones. 

The  concept  of  feedback  control  system  coupled  with 
optimization  procedure  to  design  recoil  mechanisms  was 
demonstrated  to  be  an  efficient  and  very  effective  tool. 

The  flexibility  of  feedback  control  is  added  retaining 
robustness  of  design. 

The  techniques  of  feedback  control  and  optimization 
procedures  are  recommended  for  design  of  counter  recoil 
and  other  design  problems  in  recoil  mechanism. 
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Appendix  I 

Phase-plane-delta  simulation  of  non-linear  systems 

Phase-plane-delta  method  is  an  efficient  and  fast 
simulation  method  for  non-linear  second  order  differential 
equations.  A generalized  phase  plane  method  is  also 
suitable  for  higher  order  non-linear  differential 
models;  for  details  see  "Mechanical  Systems  Analysis" 
by  A.  Seireg  [ 1 ] . The  procedure  is  discussed  below 

for  second  order  equation  reproduced  with  some  variation 
from  Prof.  Seireg' s book. 

Consider  a linear  homogeneous  second  order  differ- 
ential equation 

X + Wj^^x  = 0 (1) 

The  solution  to  this  equation  being 

X = A sin  (Wj^t  + (j))  (2) 

X = Aw^  cos  (Wj^t  + <())  (3) 


Prom  equation  2 and  3,  it  can  be  seen  that 


x2  + (i_)2 


= A" 


which  represents  a circle  in  the  x versus  ^ plane. 

This  type  of  plot  is  known  as  the  phase  plane  and  is  a 
very  powerful  and  useful  tool  in  the  analysis  of  dyna- 
mic problems.  It  could  be  seen  from  Fig.?6  that  the 
horizontal  projection  of  the  motion  on  a time  scale  will 
give  the  velocity  curve  plotted  to  a scale  w^^.  The 


Figure  26  Phase  Plane  Method 


x/p 


Figure  27  Phase  Plane  Delta  Method 


starting  point  Xq  on  the  circle  is  defined  by  initial 
conditions  at  t = 0,say  x = Xq  and  x = Xq,-  the 
amplitude  of  motion  being 

A = 

The  complete  procedure  is  detailed  step  by  step  as  follows 

1)  Calculate  the  natural  frequency  of  the  system 

2)  Plot  the  two  perpendicular  axes  representing  x 
and  (x/w^)  respectively.  These  two  axes  represent 
the  phase  plane. 

3)  The  origin  0 of  the  phase  plane  represents  the 
equilibrium  condition  of  the  system. 

4)  Knowing  the  initial  conditions  of  the  motion  in 
terms  of  initial  displacement  and  velocity,  a 
point  Xq  can  be  plotted  in  the  phase  plane. 

5)  The  free  vibration  of  the  system  is  repre- 
sented by  a circle.  The  center  of  the  circle  is 
the  equilibrium  point  0.  The  radius  of  the 
circle  is  the  distance  OXq  joining  the  origin  to 
the  initial  point  Xq. 

6)  The  conditions  of  motion  at  any  time  t^^  are 

represented  by  the  coordinates  of  the  point  X^^ 
on  the  circle  at  a radial  angle  from  Xq-  in 

the  counterclockwise  direction. 


7) 


The  projection  of  the  phase-plane  plot  on  a time 
scale  in  the  direction  of  x,  ( — ) gives  the  dis^ 
placement  and  velocity  functions  respectively. 

8)  All  information  concerning  the  time  history  of 
the  motion  can  be  obtained  fron  the  phase-plane 
plot. 

This  scheme  can  be  implemented  on  a digital  com- 
puter 

The  phase-plane-delta  method 

The  phase  plane  method  previously  discussed  is  used 
to  study  free  vibrations  of  a linear  system.  With  a 
slight  modification,  phase-plane  method  can  be  used  for 
forced  non-linear  vibration .problems  as  follows. 

Let  the  system  be 

mx  + k f(x,x,t)  - F(t)  = 0 (4) 

Rewriting  the  equation 

X + — [x  + f (x,x,t)  - X - F(t)/k  ] = 0 

m 

or 

X + p^  [x  + 6x]  = 0 (5) 

where  6x  = [f(x,x,t)  - x - F(t)/k  ] & P = /kTm 

and  is  a known  function  of  x,  x,  t,  F(t).  Equation  5 

2 

takes  the  same  form  as  equation  1 where  p takes  the 
2 

place  of  Wj^  and  x has  to  be  continually  modified  by  an 
amount  6x  which  is  a fxinction  x,  X,  t,  and  F(t). 
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Therefore,  at  any  instant  of  time  the  nonlinear  equation  can 
be  represented  in  the  phase-plane  as  a free  vibration  of 
a linear  system  with  a continually  changing  datvim. 

The  procedure  is  outlined  step  by  step,  as  follows. 

1)  Write  the  equation  in  the  form 

X + p^  (x  + 6x)  = 0 

2 

where  p is  a constant. 

2)  Express  6x  as  a function  x,  x,  t,  F(t). 

3)  Plot  the  phase  plane  axes  to  represent  x,  (x/p) 

Fig.  27. 

4)  Knowing  initial  conditions  Xq,  Xq,  a point  Xq  can 
be  plotted  in  the  phase  plane. 

5)  The  value  of  (^x)q  is  also  calculated  for  the  initial 
condition  at  ' time  t = tg. 

6)  The  center  Oq  of  the  instantaneous  linear  vibration 
at  the  initial  phase  of  the  motion  ia  therefore 
located  on  the  x-axis  at  a distance  (5x)q  from  the 
origin.  When  6 is  positive,  the  distance  is  taken 
on  the  negative  x-axis  and  vice-versa. 

7)  The  instantaneous  free  vibration  at  the  initial 
phase  of  the  motion  can  be  represented  by  a small 
arc  of  a circle  with  center  Oq  and  radius  OqXq  in 
the  counter-clockwise  direction. 

8)  A new  condition  X^^  is  reached  after  a small  increment 
of  time  tj^  where  the  displacement  x^^  and  the  velocity 
is  x^. 


9)  The  value  of  (&x) ^ corresponding  to  this  new 
condition  can  be  cSilculated,  and  a new  center 
of  oscillation  is  determined. 

10)  The  new  condition  corresponding  to  an  additional 
increment  of  time  t2  can  be  obtained  from  the  arc 
of  a circle  with  center  having  a radius  0^?! 
and  subtending  an  angle  pt2  in  the  counterclockwise 
direction. 

11)  This  procedure  is  then  repeated  until  the  time 
required  is  reached. 

12)  It  should  be  noted  here  that  the  accuracy  of  this 
procedure  depends  on  the  magnitude  of  the  angular 
displacement  pt.  Better  accuracy  can  be  attained 
also  by  iteration.  This  means  that  after  is 
determined,  a value  of  6x  is  calculated  which 
corresponds  to  a point  midway  between  X^  and  X^. 
Using  this  new  value  of  6x,  a better  approximation 
for  the  new  position  after  a small  increment  of 
time  can  be  obtained.  This  procedure  can  be 
repeated  until  a desired  accuracy  is  achieved. 

13)  This  procedure  is  suitable  for  programming  on  a 
high-speed  digital  computer. 
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APPENDIX  II 

ALGORITHMS  FOR  OPTIMIZATION 
OF  NON-LINEAR  FUNCTIONS 

Three  algorithms  for  non-linear  optimization  are  discussed: 

1)  Davidon-Fletcher-Powell  unconstrained  optimization  with 
cubic  interpolation. 

2)  Davidon-Fletcher-Powell  unconstrained  optimization  with 
linear  search  by  golden  section. 

3)  Constrained  optimization  by  Random  search. 

A short  introduction  to  gradient  methods  of  optimization  is  presented. 
The  problem  of  non-linear  optimization  can  be  stated  as 
Find  a vector  x of  parameters  to  minimize 

J(x),  an  objective  function  subject  to  certain  con- 
straint of  type 

C^.(x)  ^ , = or  ^ b. 

We  can  express  J(x)  in  Taylor  series  expansion  as 

J(x+  Ax)  = J(x)  + Ax'g  + j Ax'  GAx  + (1) 

where  Ax  is  a vector  of  increment  in  x 
g is  the  gradient  vector 

G is  the  Hessian  Matrix  of  second  order  partial 
derivatives 

The  first  order  approximation  is 


J(x+Ax)  = J(x)  + Ax'g 
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The  reduction  in  the  function  J(x)  for  moving  to  x+Ax  is  Ax' g and  is 
maximum  if  we  move  in  the  negative  direction  of  gradient  vector  g. 


Ax  = “X  -2- 

|g| 


(2) 


The  optimum  value  of  X is  found  by  univariate  linear  search  in  the 
direction  — ^ . 

|g| 

This  is  called  the  steepest  descent  method.  This  method  converges 
very  slowly.  The  extension  of  this  method  is  the  conjugate  gradient 
method.  Assuming  the  objective  function  is  quadratic  of  the  form. 


J(x)  = y x’Gx  + bx  + c 


(3) 


one  can  find  a sequence  of  search  directions  u\  u ,...,u.  such  that 

12  K 


u,  Guj  = 0 


l7j 


(4) 


By  moving  sequentially  in  the  directions  Uj.Uj,. 
variate  searches  we  can  reach  the  optimum 


X*  . = X + e X.  u. 

im  n 0 11 

i=l 


,U|^  with  linear  uni- 


(5) 


The  steepest  descent  and  conjugate  gradients  are  illustrated  in 
Fig.  28  for  two  parameter  case 
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Steepest  Descent 


Conjugate  Gradient 


Figure  28a 


Figure  28b  Concept  of  Conjugate 
Gradient 
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The  conjugate  gradient  method  transforms  concentric  elliptic 
contours  into  concentric  circular  contours.  The  circular  contour  has 
the  desirable  property  of  the  normal  to  the  tangent  at  any  point  passes 
through  the  center.  Once  the  direction  of  normal  is  found  in  this  plane, 
it  can  be  transformed  back  to  the  original  plane  and  direction  of  the 
increment  in  x can  be  found  out  and  this  will  pass  through  the  minimum 
for  the  quadratic  function. 

(3)  Newton-Raphson  Method 

J(x+Ax)  = J(x)  + Ax'g  + j Ax'G  Ax  (6) 

differentiating  with  respect  to  Ax 

. g . GA. 

for  x+Ax  to  be  a optimum  point 

3J(x+Ax)  _ rt 
3Ax  ^ 

g + GAx  =0  , Ax  = -G~^g  (7) 

This  method  is  called  Newton-Raphson  and  is  a very  efficient  method 
if  6 and  g are  available. 


The  problem  for  non-quadratic  function  minimization  is  to  find  G 
which  is  positive  definite  and  the  question  of  how  good  an  approximation 
of  second  order  Taylor  expansion  is.  The  combination  of  conjugate 
gradient  method  for  fast  movement  when  away  from  minimum  can  be  utilized 
with  the  conjagacy  being  with  respect  to  the  Hessian  G.  The 
Davidon-Fletcher-Powel 1 (DFP)  algorithm  does  exactly  that. 
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The  DFP  algorithms  continuously  updates  H,  the  inverse  of  the  hessian 
matrix  with  linear  search  in  direction  of  conjugate  gradient  by  cubic 
interpolation.  The  step  by  step  procedure  given  below  is  reproduced 
with  variation  from  [4].  The  details  can  be  found  in  [2,3,4  & 5]. 

1)  Set  = I and  let  k be  the  current  iteration  number,  then  set 

dk  = -Hk9k  (8) 


then,  d|^  is  the  direction  of  search  from  the  current  point 

2)  Perform  a linear  search  to  find  X |^(>0),  where  is  the  value 
of  X|^  that  minimizes  ^ 

3)  Set  Ax,^  = X*^d^  (9) 

4)  Set  + AX|^  (10) 

giving  the  new  current  point 

5)  Evaluate  and  gj^^^ 

6)  Set  Ag,^  = g|^^^  - g,^  (11) 


7)  H 


k+1  ^k 


^k  ^'k  , ¥9k^g’k^k 


Ax'kAOk  A'g^H^Ag^ 
8)  Set  k=k+l  and  return  to  step  1 . 


9)  Stop  when  either  |d|^|  or  every  component  of  dj^  is  smaller  than 
some  prescribed  amount. 

The  linear  search  in  step  2 can  be  performed  in  two  ways  with  cubic 
interpolation  or  golden  section  search. 
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Cubic  interpolation  procedure; 

1)  Evaluate  = J(xj^)  and  6^  = dj^ 


(12) 


check  < 0.  Compute  a by 


a = min[2, g ] 


2(Jo-Jg) 


(13) 


0 


★ 

where  Je  is  the  estimated  values  of  J(xj^+X 


2)  Evaluate  = J(xj^+a  dj^)  and  = g'dj^ 

3)  If  6^  >0  or  If  > Jq’  P>"oceed  to  rule  5 otherwise  go 
to  rule  4. 

4)  Replace  a by  2a,  return  to  rule  2. 

5)  Interpolate  in  the  interval  [0,a]  for  using 


6)  Return  to  rule  5 to  repeat  the  interpolation  in  the  smaller 
interval  [0,X*j^]  or  [Xj^,a]  accordingly  as 


stop  when  the  interval  of  interpolation  has  decreased  to 
some  prescribed  value. 


a ~ 6„-6.+2w 
a 0 


(14) 


where 


1/2 


w = (z^-GjGj) 


(15) 


and 


9'm  <*!(><><>'•<» 


J. 
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The  complex  method  of  M.J.  Box  [6],  [7],  with  random  search  generates 
a simplex  of  n+1  points  where  n is  the  number  of  parameters  to  be  opti- 
mized. The  function  is  computed  at  each  of  the  vertexes  and  the  worst 
vertex  is  discarded.  The  reflection  of  this  discarded  vertex  into  the 
centroid  of  the  remaining  vertices  is  taken  as  a new  vertex  and  search 
is  carried  out  in  two  dimensions,  the  nrocedure  is  as  follows: 


Three  initial  vertices  are  x^,  x^,  Let  x^  be  the  worst  vertex 
and  x^  be  the  centroid  of  Xj  and  x^.  Then  x,  is  computed  by  the  rela- 
tion 


where  is  vertex  with 
repeated  until  minimum  is 


minimum  function  value, 
obtained. 


The  procedure  is 
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Appendix  III 

Design  Data  for  M-.37  Recoil  Mechanism 

M-37-  recoil  mechanism  is  designed  for  firing  zone 

7.  All  the  pertinent  parameters  necessary  for  simulation 
of  recoil  mechanism  model  developed  in  Chapter  1 are 
given  below.  The  variable  area  of  the  groove  machined  in 
the  floating  piston  is  tabulated  in  Table  2.  Breech  forces 
are  tabulated  in  Tables  3 through  7 for  zones  1,  5 through 

8.  These  breech  forces  are  simulated  breech  force  data 
rather  than  actual  test  data.  All  the  design  data  listed 
here  was  provided  by  Control  and  Stabilization  group  of 
General  Rodman  Laboratory,  Army  Weapons  Command  at  Rock 
Island,  Illinois. 

2 

irij^  - Mass  of  recoiling  parts  = 3.6658  Ibf  sec  /in 

2 

nip  - Mass  of  floating  piston  = .06735  Ibf  sec  /in 
p - Mass  density  of  hydraulic  fluid  = .78E-4 
Pq  - Initial  gas  pressure  = 1153  psi 

3 

Vq  - Initial  volume  of  gas  = 513  in 
R - Gas  constant  = 1.68 

2 

Ap  - Recoil  piston  area  = 2.9906  in 

2 

A^  - Control  rod  area  = 2.4053  in 

2 

Ajj  - Floating  piston  area  = 11.781  in 

Cj^=C2*=C2  - Discharge  coefficients  = . 8 

2 

Aj^  - Area  of  orifice  1 = 1.25  in 

2 

A2  “ Area  of  orifice  2 = .4536  in 
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TABLE  2- 


AREA  nF  VARIABLE  HRIFIcE  A3 


X 

A3(X) 

X 

A3(X) 

0.00000 

.08680 

15.12600 

.07900 

.50A20 

.09370 

15.63020 

.07710 

1.00840 

.09410 

16.13440 

.07520 

1.5X260 

.09630 

16.63860 

.07320 

2.01680 

.09950 

17.14280 

.07130 

2.52100 

.10110 

17.64700 

.06940 

3.02520 

.10210 

18.15120 

.06750 

3.52940 

.10250 

18.65540 

.06540 

4.03360 

.10270 

19.15960 

.06320 

4.53780 

.10290 

19.66380 

.06100 

5.04200 

.10280 

20.16800 

.05880 

5.54620 

.10250 

20.67220 

.05630 

6.05040 

.10190 

21.17640 

.05380 

6.55460 

.10120 

21.68060 

.05120 

7.05880 

.10060 

22.18480 

.04870 

7.56300 

.09960 

22.68900 

.04580 

8.06720 

.09870 

23.19320 

.04290 

8.57140 

.09770 

23.69740 

.03980 

9 . 07560 

.09670 

24.20160 

.02650 

9.57980 

.09540 

24.70580 

.02270 

10.08400 

.09410 

25.21000 

.02870 

10.58820 

.09290 

25.71420 

.02410 

11.09240 

.09150 

26.21840 

.01910 

11.59660 

.09000 

26.72260 

.01330 

12.10080 

.08850 

27.22680 

.0C620 

12.60500 

.08720 

27.73100 

O.OCOOO 

13.10920 

.08570 

28.23520 

O.OCOOO 

13.61340 

.08410 

28.73940 

O.OCOOO 

14.11760 

.08250 

29.24360 

O.OCOOO 

14.62180 

.08090 

J 
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table  3-  r.REECM  PHRcE  FOR  zn^E  I 


T 

8(T) 

T 

B(TJ 

OtOOOOO 

0.0 

,03800 

1757,0 

,00050 

5.0 

,04000 

1526,0 

,0010:'. 

6.0 

,04200 

1329,0 

,00150 

7,0 

,04400 

1159,0 

,00200 

8,0 

,04600 

1013,0 

.00250 

0.0 

,04000 

887,0 

,0030« 

10,0 

.05000 

778,0 

,00350 

II. 0- 

.05200 

683,0 

,00400 

12.0 

.05400 

6C2,0 

,00450 

10.0 

,05600 

530,0 

,00500 

72.0 

,05800 

468,0 

.0055* 

1177. 0 

.06000 

414,0 

.00600 

2041.0 

.06200 

367,0 

.0065". 

3287.0 

,06400 

326,0 

,00700 

5272.0 

,06600 

289.0 

.00750 

3393.0 

.06000 

257,0 

,00800 

13220.0 

,07000 

229,0 

.00050 

20477.0 

,07200 

2C5.0 

.onooo 

30024.0 

.07400 

183,0 

,00957 

.44990.0 

,07600 

163,0 

,01000 

62180.0 

,07800 

146.0 

.oioso 

30403.0 

,08000 

131,0 

.oiloo 

96098.0 

,06200 

118,0 

.0115'. 

105714,0 

,06400 

1C6.0 

.Oils*' 

l«75nn,9 

,08600 

95,0 

,0125:i 

102620,0 

,08800 

86,0 

.01300 

93260,0 

,09000 

77,0 

.Ol35r 

82047,0 

.09200 

70,0 

.01400 

70792.0 

,09400 

63,0 

.01450 

60482.0 

.09600 

57.0 

.01500 

51500,0 

,09800 

52,0 

.01550 

43889,0 

,10000 

47,0 

,01600 

37532.0 

,10200 

43.0 

,01651 

32252,0 

,10400 

39,0 

.01700 

27372,0 

,10600 

35,0 

.01750 

24231.0 

.10800 

32.0 

.OlfiOO 

21192.0 

,11000 

29,0 

.01850 

13643.0 

tll200 

27,0 

.01900 

16493,0 

.11400 

24,0 

.01950 

14669,0 

•11600 

22,0 

.02000 

13113,0 

,11800 

20,0 

,02050 

11778,0 

. 12000 

19,0 

.02100 

10626,0 

.12200 

17,0 

.02150 

9626.0 

,12400 

16,0 

.02200 

8755.0 

.12600 

14.0 

.0225*' 

7992,0 

.12800 

13,0 

,02300 

7320.0 

.13000 

12,0 

,02350 

9726,0 

.13200 

U.O 

,02400 

6199.0 

.13400 

10,0 

,02450 

5729.0 

,13600 

9.0 

,02500 

5309,0 

•13800 

9.0 

.02550 

4933.0 

. 14000 

8.0 

.02600 

4593.0 

.14200 

7.0 

.02650 

4286,0 

, 14400 

* f w 

7.0 

.02700 

4009,0 

, 14600 

• t w 
6*0 

.02724 

3884,0 

,14800 

6.0 

.02750 

3008.0 

. 15000 

5,0 

.0280.0 

3665,0 

.15200 

5,0 

,03000 

3149.0 

.20000 

0.0 

.03200 

2712.0 

,03400 

2342.0 

,03600 

2026,0 
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TABLE  4-  RREECH  FURCE  FHR  ZHNC  5 


T 

8(T) 

T 

B(T) 

0.00000 

159.0 

,03600 

2426.0 

.00050 

246,0 

,03800 

2069,0 

.00100 

513,0 

,04000 

1769,0 

.00150 

938,0 

,04200 

1517.0 

.00200 

1724,0 

,04400 

13C4.0 

.00250 

3118,0 

,04600 

1123,0 

.00300 

5642,0 

,04800 

969.0 

.00350 

10127,0 

,05000 

839,0 

.00400 

18001.0 

,05200 

727,0 

.00450 

31434,0 

,05400 

632,0 

.00500 

53257.0 

,05600 

550,0 

.00500 

85873,0 

,05800 

480,0 

.00600 

128386.0 

,06000 

419,0 

.00650 

173036,0 

,06200 

367,0 

.00700 

206142.0 

,06400 

322.0 

.0074 A 

217035,0 

,06600 

283,0 

.00750 

217008,0 

,06800 

249,0 

,00800 

205962,0 

,07000 

220,0 

.00850 

181687,0 

,07200 

194.0 

,00900 

153385,0 

,07400 

172.0 

,00950 

126693.0 

,07600 

152,0 

.01000 

103872,0 

,07800 

135.0 

,01050 

05258.0 

,08000 

120,0 

,01100 

70386.0 

,08200 

IC7.0 

.01200 

49196,0 

.08400 

95,0 

,01250 

41695.0 

,08600 

85,0 

.01300 

35654,0 

,08800 

76,0 

.01350 

30747,0 

.09000 

68.0 

.01400 

26725,0 

,09200 

61.0 

,01450 

23400,0 

,09400 

54.0 

.01500 

20628,0 

.09600 

49,0 

.01550 

18298,0 

.09800 

44,0 

.01600 

16324.0 

,10000 

39,0 

.01650 

14641,0 

,10200 

36,0 

.01700 

13195.0 

,10400 

32.0 

,01734 

12315,0 

.10600 

29,0 

.01750 

12135.0 

,10800 

26.0 

.01800 

11575,0 

.11000 

24,0 

,02000 

9605,0 

,11200 

21.0 

.02200 

7998,0 

,11400 

19.0 

.02400 

6683,0 

.11600 

18.0 

,02600 

5602.0 

.11800 

16,0 

.02800 

4711,0 

,20000 

0,0 

,03000 

3973.0 

.03200 

3361,0 

,03400 

2852.0 

table  5-  BREECH  FllRCE  FPR  2PNE  6 


T 

R(T) 

T 

E (T) 

0.00000 

384.0 

.03500 

2545,0 

.00050 

706,0 

,03700 

2135,0 

.00100 

1491,0 

,03900 

1797,0 

.00150 

2937,0 

.04100 

1517.0 

.00200 

5000,0 

.04300 

1285,0 

.00250 

11343,0 

.04500 

1091,0 

.00300 

21947,0 

.04700 

929,0 

.00350 

41614,0 

.04900 

793.0 

.00400 

76017,0 

.05100 

678,0 

.00450 

130033.0 

.05300 

582,0 

.00500 

200301,0 

.05500 

5C0.0 

.00550 

267020,0 

.05700 

431,0 

.00600 

302300,0 

.05900 

373,0 

.00616 

304395,0 

.06100 

323,0 

.00650 

295349,0 

.06300 

280.0 

.00700 

259560,0 

.06500 

243,0 

.00750 

214473,0 

.06700 

212.0 

.00800 

172277,0 

.06900 

185.0 

.00050 

137320.0 

.07100 

162,0 

.00900 

109863,0 

.07300 

142.0 

.00950 

88691,0 

.07500 

124.0 

.01000 

72415.0 

.07700 

1C9.0 

.01050 

59035.0 

.07900 

96,0 

.01100 

50018,0 

.08100 

85.0 

.01150 

42274.0 

,08300 

75.0 

.01400 

2088t3.0 

,08500 

66,0 

.01431 

19380,0 

,08700 

59.0 

.01450 

18986,0 

,08900 

52.0 

.01500 

17983,0 

,09100 

46,0 

.01700 

14514.0 

.09300 

41.0 

.01900 

11769,0 

.09500 

37,0 

.02100 

9586,0 

.09700 

33,0 

.02300 

7841.0 

.09900 

29,0 

.02500 

6440,0 

,10100 

26,0 

.02700 

5310.0 

.10300 

24,0 

.02900 

4395.0 

,20000 

0.0 

.03100 

3651,0 

.03300 

3043,0 
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TAI3LE  6-  BRECCil  HlKCC  FOP.  ZH  IE  7 


T 

B(T) 

T 

3(T) 

o.noooo 

2809,6 

,01700 

13971.0 

.00020 

4589.6 

,01740 

13354.0 

.00060 

9058,6 

,01780 

12765,0 

.00100 

18069,0 

.01820 

122C4.0 

.00140 

35565,0 

,01860 

11668.0 

.00180 

68543,0 

.01900 

11156,0 

.00220 

126940.0 

,01940 

10674,0 

.00260 

219070.0 

.01980 

10214.0 

.00300 

337720.0 

,02020 

9775,1 

.00340 

446170,0 

,02060 

9355.5 

.00380 

496900,0 

,02X00 

8954.5 

.00420 

476070,0 

,02140 

8576,1 

,00460 

411400.0 

.02180 

8215.1 

.00500 

335660,0 

.02220 

7869,6 

,00540 

267300,0 

.02226 

7538,8 

.00580 

211870,0 

,02300 

7222.5 

.00620 

168900,0 

.02340 

6924.1 

.00660 

136050,0 

,02380 

6638.1 

.00700 

110950,0 

,02420 

6364,5 

.00740 

91613,0 

,02460 

61C2.9 

,00780 

76562.0 

.02500 

5852.1 

,00820 

64706,0 

,02540 

5615.3 

,00860 

55253.0 

.02580 

5389,0 

.00900 

47627.0 

,02620 

5171.7 

,00940 

41322.0 

,02660 

4963.0 

.00980 

36281,0 

.02700 

4763.3 

.01020 

32030,0 

,02740 

4574,0 

,01060 

29704.0 

,02700 

4393.7 

.01100 

28358.0 

.02820 

4220.4 

,01140 

27013.0 

,02860 

4054.0 

.01180 

25736,0 

,02900 

3894.2 

.01220 

24521.0 

,02940 

3743,0 

.01260 

23366,0 

.02900 

3598.0 

.01300 

22263.0 

.03020 

3458,6 

.01340 

21237.0 

.04020 

1361,8 

.01380 

20257.0 

.05020, 

578.4 

.01420 

19324,0 

,06020 

263.3 

,01460 

18435.0 

.07020 

126.1 

,01500 

17588,0 

,08020 

63.5 

.01540 

16794.0 

,09020 

0,0 

.01580 

16036,0 

.2000C 

0,0 

.01620 

15314,0 

.01660 

14626,0 
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table?-  breech  force  for  Z0‘1E  8 


T B(T)  T B(T) 


OiOOOOO 

3090,6 

.0002-'' 

5048,6 

.00060 

'»964,5 

.00100 

19075,9 

.00140 

39X21.5 

.00180 

75397,3 

.00220 

139634.0 

.00260 

240977,0 

.00300 

371491.9 

.00340 

490786.9 

.00380 

546589.9 

.00420 

523676,9 

.00460 

452539,9 

.00500 

369225,9 

.00540 

294029,9 

.0058<^ 

233057,0 

.00620 

135790,0 

,00660 

149655,0 

,00700 

122045.0 

.00740 

100774.3 

.00780 

84218,2 

.00820 

71176.6 

.00860 

60778,3 

.00900 

52389.7 

.00940 

45454.2 

.00980 

39909,1 

.0X020 

35233.0 

.0X060 

32762.4 

.0X100 

31193,8 

.01140 

29714.3 

.0X180 

28309,6 

.0X220 

26973.1 

.0X260 

25702.6 

.0X300 

24494,8 

.0X340 

23360,7 

.0X380 

22282,7 

.0X420 

21256,4 

.0X460 

20278.5 

.0X500 

19346.8 

.0X540 

18473,4 

,0X580 

17639,6 

.0X620 

16845.4 

.0X660 

16088.6 

.01700 

15368.1 

.01740 

14669.4 

,01730 

14041.5 

.01820 

13424.4 

.01860 

12834,6 

,01900 

12271.6 

.01940 

11741.4 

.01980 

11235.4 

.02020 

10752.6 

,02060 

10291.0 

,02100 

9849.9 

.02140 

9433.7 

,02180 

9036.6 

,02220 

8656,6 

,02226 

8292.7 

,02300 

7944.7 

,02340 

7616,5 

,02380 

73C1.9 

,02420 

70C0.9 

,02460 

67X3,2 

.02500 

6437.3 

.02540 

6176.8 

.02580 

5927,9 

.02620 

5668,9 

,02660 

5459.3 

.02700 

5239.6 

.02740 

5031,4 

.02780 

4833.1 

.02820 

4642,4 

,02860 

4459.4 

,02900 

4263,6 

,02940 

4117.3 

,02980 

3957,8 

,03020 

38C4.5 

,0402c 

X498.0 

,05020 

636,2 

,06020 

289,6 

,07020 

140.9 

,00020 

69,8 

,09020 

0,0 

,20000 

0,0 
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